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ABSTRACT 

Density inhomogeneities along the line-of-sight distort fluctuations in the cosmic 
microwave background. Usually, this effect is thought of as a small second-order 
effect that mildly alters the statistics of the microwave background fluctuations. We 
show that there is a first-order effect that is potentially observable if we combine 
microwave background maps with large redshift surveys. We introduce a new 
quantity that measures this lensing effect, < T[50 • VT) >, where T is the microwave 
background temperature and 59 is the lensing due to matter in the region probed by 
the redshift survey. We show that the expected signal is first order in the gravitational 
lensing bending angle, < (66)'^ >^/^, and find that it should be easily detectable, 
(S/N)~ 15 — 35, if we combine the Microwave Anisotropy Probe satellite and Sloan 
Digital Sky Survey data. Measurements of this cross-correlation will directly probe the 
"bias" factor, the relationship between fiuctuations in mass and fiuctuations in galaxy 
counts. 
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1. Introduction 



The next several years should be very exciting for cosmologists: Microwave Anisotropy 
Probe (MAP; [Wright et al. 1996| ) and PLANCK ( pouchet et al. 1995D will make high resolution 
maps of the microwave background sky; while the Sloan Digital Sky Survey (SDSS; |Gunn and 
Weinberg 1995 ; see also http:/ /www- sdss.fnal.gov:8000/| ) will measure redshifts of 10^ galaxies 
and positions of 10*^ galaxies. In this paper, we explore the direct connection between these two 
measurements through gravitational lensing: the path of a cosmic microwave background (CMB) 
photon is distorted by inhomogeneities in the matter distribution; galaxy surveys detect these 
inhomogeneities as fluctuations in galaxy number counts. 

The effect of the gravitational lensing on the CMB anisotropies has been studied by many 
authors. The uncomfortably low upper limits (Uson and Wilkinson 1984, Readhead et al. 1989| ) 
provoked a great deal of controversy (Kashlinsky 198J, Tomita 1988, Sasaki 198£, Watanabe and] 
p?omita 1991 ) about the possibility that gravitational lensing washes out the intrinsic fluctuation. 
After the detection by the Cosmic Background Explorer (COBE; pmoot et al. 199^ ), there has 
been renewal of interest ( Linder 1990a , b, Cayon et al. 1993a| , b, [Seljak 1996| ) in investigating 
how the CMB power spectrum is redistributed owing to gravitational lensing. For example, 
Seljak (1996) recently presented detailed calculations of gravitationally deflected CMB power 
spectra, including the effect of the nonlinear evolution of matter density fluctuations. His result 
shows, however, that the modification of the CMB power spectrum is a second-order effect of the 
photon bending angle and less than a few percent on angular scales greater than ten arcminutes. 
Hence, the lensing effect on the CMB spectrum itself is extremely difficult to detect, even with 
observations such as the MAP project. Linder (1997) has also studied the effects of lensing on 
the correlation function and has introduced a cross-correlation function similar to the one that we 
study here. 

In this paper, we introduce a cross-correlation function that is sensitive to the gravitational 
lensing correlations between the temperature fluctuations and matter density fluctuations. We 
show that the cross-correlation is first-order in the bending angle so it should be easier to detect if 
we have both accurate CMB maps and redshift surveys. We quantitatively estimate its magnitude 
and its cosmic variance in cold dark matter (CDM) universes. The rest of the paper is organized 
as follows. We review the formalism developed by Seljak (1996) for computing the angular 
excursion of the CMB photon paths on celestial sphere in section 2. In section 3, we formulate 
the cross-correlation between matter density inhomogeneities and CMB temperature fluctuations. 
Section 4 concludes. 



2. Gravitational Lensing 



In this section, we review gravitational lensing by density fluctuations. We follow the power 
spectrum approach of Seljak (1994, 1996). We focus on the angular excursions produced by matter 
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fluctuations at low redshifts, where they can be most easily inferred from redshift surveys. 

Fluctuations in matter density, S, generate variations in the gravitational potential, 

V2(/> = A7rGpba^5, (1) 

where G is the gravitational constant, and ph is the mean background mass density. Conventionally, 
the matter density fluctuations are related to the fluctuations in galaxy counts by a linear biasing 
parameter, b: 

6g=b5. (2) 

Since most of the lensing effects will be produced by fluctuations on large physical scales 
{k < 0.1/i~^Mpc), the linear biasing model will hopefully be valid. It is important to note that 
detailed nonlinear and/or time-dependent biasing may somewhat change the statistics we present 
in this paper. 

A photon emitted at some angular position 6 has been deflected by gravitational lensing 
during its long travel, with the result that it is observed at different angular position, tp. The 
photon angular excursion on celestial sphere is given by Seljak (1994): 

- V = = -2 r^'^ dx' W{x', Xdec) Vxc/', (3) 

Jo 

where is transverse component of the potential gradient with respect to the photon path, 

W{x, Xdcc ) = sin^^ {xdec - X) / sinK (Xdec ) (4) 

is a projection operator on celestial sphere, and Xdec = x(^doc) is unperturbed comoving radial 
distance corresponding to redshift Zdec at decoupling time. In equation (Q), smxiu) = sm{u),u, 
and sinh(u) in a closed, flat, and open universe, respectively. 

Next, we consider the relative angular excursion 56 — 56' of a photon pair emitted from 
angular positions 6 and 6' . We restrict our calculation to the small angular separation limit, 
^ = \6 — 6'\ <^ \, and assume that the relative angular excursion 56 — 56' obeys Gaussian statistics. 
Lensing is primarily due to scattering events from mass fluctuations on the 10 - 100 Mpc scale. As 
there are 30 - 300 of these fluctuations between the surface of last scatter and the present along 
each photon path, the central limit theorem implies that this is a good approximation. Following 
Seljak (1994), we characterize the statistics of the lensing fluctuations by its root-mean-square 
dispersion: 

a{i-z) = 2-V2([50(z)-50'(z)]^)^^' = [Ql(O;z)-Cgl(C;z)]^/^ (5) 
Qi(C;^) = -/ k^dk dx'P4k,To-x')W\x',Xdcc)Mk^smKx'), 

TT Jo Jo 

where ( )^ denotes the averaging over pairs observed with fixed angular separation ^, Jq is the 
Bessel function of order 0, and P^{k) is the gravitational potential power spectrum. The power 
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Fig. 1. — The dispersion of the relative angular excursion, a{6;z). la: Dependence on the 
separation angle 9. Solid lines show the standard cold dark matter (SCDM) model with = 1; 
^vo = 0, h = 0.5, and cg = 1.2. Broken lines show the cosmological constant dominated cold dark 
matter (ACDM) model with 0,^0 = 0.3, O^^o = 0.7, h = 0.7, and as = 1-0. For each model, the 
curves show the case in which z = ^deo Ij 0.5, and 0.2, from upper to lower, lb: Dependence on 
the redshift z at 9 = 0.21 degree in the SCDM model. 



spectrum of the potential fluctuations are related to the power spectrum of the density fluctuations 
through, 

P^k, r) = {9/A)^l{T)H\T)a\T)k-*P{K r), (6) 

where Q.m is the mass density parameter given by = 8TTGpb/{3H'^), and H{t) is the Hubble 
parameter. 

Figures la shows a{9; z) as functions of 9 for several redshift values. We consider throughout 
this paper two cosmological models: one is the standard CDM (SCDM) model with ^7^0 = 1; 
^vO = 0, h = 0.5, and erg = 1.2, and the other is a low-density, cosmological constant dominated 
CDM (ACDM) model with n^o = 0.3, n^o = 0.7, h = 0.7, and erg = 1.0, the best fit model of 
Ostriker and Steinhardt (1995). Here, ^ImO ^-nd fi^o are the present mass density and the present 
vacuum energy density normalized by the critical density; h is the present Hubble parameter in 
units of 100kms~^Mpc~^; ug is the mass fluctuations within a sphere of radius 8 Mpc. We use 
the COBE normalized value (Bunn and White 1997) for ug. In numerical calculation of a(^;z), 



we have used the fitting formula for CDM linear transfer function given in Bardeen et al. (1986). 

The SDSS project measures the redshift of galaxies up to z ~ 0.2 within solid angle 
^SDSS = vrsteradian. Furthermore, we expect to obtain photometric redshifts (Connolly et al. 



1995]) of galaxies up to z ~ 1 with fairly small uncertainties. Then we will obtain the galaxy 
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number density perturbation, Sg{r), within solid angle fisDSS for z < 0.2, and for ^; < 1 with 

some uncertainties due to error in photometric redshifts. We can see from Figures lb that the 
matter within z < 0.2, z < 0.5, and z < 1 contributes 15%, 35%, and 55%, respectively, to the 
angular excursion from the last scattering surface. Hereafter, the values with subscript A denote 
the integral from the observer to some redshift z, the edge of the survey data, and those with 
subscript B denote the integral over remaining part; for example: 

66 A = -2 r^'\x'W{x',Xdec)'V±<P, (7) 
Jo 

fXdec 

50 B = -2 / dx'W(x',Xdec)Vi<^. (8) 

The total lensing deviation, 56 = 56a + 56b- Its variance is the sum of the contribution from the 
two regions: 

{{56f) = {(SOa?) + {{56b)') ^ 2{a\ + 4). (9) 

For our purposes, the distant lensing is unimportant; its only effect is to slightly reduce the 
amplitude of the temperature fluctuations. 



3. Lensing the Microwave Background 

Gravitational lensing distorts the microwave background sky: 

f{il}) = T [6{il})\ =T{^ + 56)= T{iP) + 56 • VT{iP) + ^[56 • Vf T + ■ ■ ■ . (10) 

Here T denotes the measured temperature map and T denotes the unlensed temperature map. 

This distortion alters the statistics of the microwave background by smearing out temperature 
correlations (Seljak 1996): 



[ni,)T{^')) = mm')) + {[56- vr(v)] [56' -vn^')] > 

+ {T{^){56'-VfT{^'))- (11) 
This can be alternatively written in terms of the correlation function: 



Clensed(0 — C'uniensed(0 + ^ (1^) 



where ^ = — %l^' \ . 



If we have a redshift survey, then the effects of gravitational lensing of the microwave 
background can be observed more easily. Here, we introduce a new quantity, which measures 
the cross-correlation between the temperature map and the predicted lensing: 

H (0 =M{f{^) {[56a - 56'a\ -Vf (VO))^ , (13) 
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where 50a{'4') is determined from the redshift survey, and is a normahzation factor given by 



^ a^(^fwhm) It')''" {{^Tf)''' (14) 

^ ' ffwhm ^ ' ffwhm 



at a typical value of angular separation 9 = ^fwhm- 

We can evaluate this by expanding the temperature map in a Fourier series. Since the lensing 
bending angles are small, we simplify the equations by making the plane parallel approximation 
and expanding out the unlensed temperature map: 

r(0) = ^a;exp(-iZ -0), (15) 
I 

where denotes multipole moments. We can now rewrite the effects of lensing in the multipole 
expansion: 

r(^) = T(V +56) (16) 
= aj exp [— iZ-(?/> + 66)] , 



and _ 

dT{il^') 



dip' 



il'al, exp + 66')] . (17) 

l 



We can now evaluate the lensing statistic: 



H{^) = A/'^^(^a;a*,exp(-iZ-56> + iZ'.(56'') 
I I' 

X [il ■ (66 A - 66'^)] exp {-il • -0 + il'-i^') ) (18) 



Z 

X (il • {66 A - 66'a) exp [-il ■ {66a - 66'^)]) 
x{exp[-il- {66b-66'b)]). 

Since the gravitational deflections are the sum of many small scattering due to superclusters 
and voids along the line-of-sight, we can treat x = I ■ {66 — 66') as a Gaussian random variable. 
The dispersion of x is {x'^) = I'^a'^, where we have kept only the main isotropic term. The 
anisotropic term makes a subdominant contribution to the gravitational lensing effect on two-point 
auto-correlation function of CMB ( peljak 1996| , Martinez-Gonzalez et al. 1997). Though the 
anisotropic term is not too small in this statistics, the more rigorous analysis will be given in a 
subsequent paper ( ^uginohara, Suginohara and Spergel 1997| ). Then 

{ixexp{-ix)) = -J^SLJ^ L (19) 

n=0 
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E7i-ATT(2n + l)n(x2^ 



n=0 



(2n + 1) 



;2 2 

/ a exp 



Combining equations (ITq) and (|l9| 



H{0 = Mj2{aj)l^al{0exp 
I 



exp [—il • {tp — tp')] 



X exp 



(20) 



The final term represents the smearing of the microwave fluctuations due to density fluctuations 
beyond the edge of the redshift survey. This effect is a small correction term, which we will ignore 
for the remainder of the paper. 



We average equation (pO) over angle and rewrite it on the celestial sphere, then 



exp 



Pz(cose), 



(21) 



where ^ is the angular separation on the sky, Ci is the usual multipole moment, Wi = exp(— /^cr^gj^j^) 
is the window function, and (Tbeam is the beam of the detector ( Knox 1995| ). 



The cosmic variance in the cross-correlation statistic can also be estimated by taking the 
leading order term in the expansion: 



TitP) ([56 - 60'] •VT('0') 



(22) 



1 /2 

Note that the signal-to-cosmic variance scales as (H) / {H"^) , which is proportional to cta- As 
we claimed earlier, this is a first-order effect. 

We can estimate the signal-to-cosmic variance ratio by computing the predicted signal per 
multipole: 



dtpdtp'Tiitp) (deitp) - 66' [iP')) -VTi {tP ') 

dipdip'ai exp [-il ■ {ip + 66)] [il ■ {66{ip) - 66' {jP'))] 
xa*i ex.-p[il-{ip' + 66')] 



(23) 



CiWifafexp{-t'af/2), 



2_2 
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where cj; = (T^(7r//). If we include detector noise, then 

(t^) = CiWi + w-\ (24) 



where if ^ is the measure of the detector noise ( Knox 1995| ). MAP's highest frequency 



channel has a full width at half maximum of 0.21 degree and a system noise of 

= (10.5/iK)^ degree'^ = (0.18/iK)-^ steradian. By combining the three highest frequency 
channel, the MAP system noise drop to (O.ll/uK)^ steradian. This implies that the noise plus 
cosmic variance in Hi is 

where the factor of 2Z + 1 comes from averaging over all multipoles with the same I. Note that 
there is no factor of 2 in the numerator as the variance is proportional to the product of two 
uncorrelated fields, (T^ (Vr)^V rather than the more familiar (T^) - (T^)' = 2 (T^)' . bumming 



over all multipoles yields the expected signal in the microwave maps: 



-Y - (26) 

Pafexp{-Paf) 



^(2/ + l) 



l + iwQWi)-'? 



Note that an additional cross-correlation function that can be computed from the microwave 
background fluctuations and the lensing maps: 

GiO = {f{iP) [[5eAW-60'A^')] -nfi^P')))^, (27) 

where n = |^ ^'j • 'The statistics has the simplest form among possible inclusions of gravitational 
lensing angular excursion. The cross-correlation function is reduced to 

GiO=Moa\^. (28) 

However we confirmed that the expected signal-to-noise ratio, S/N, in G is lower than in H. This 
is because, unlike in H, the cosmic variance in G contains a term proportional to (r^(i/7)T^('0')). 



4. Results and Discussion 

We have computed the expected cross-correlation between the temperature fluctuations and 
the lensing bending angle (equation (^TI)) for the standard CDM model and the "best fit" vacuum 
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dominated model ( pstriker and Steinhardt 1995| ). In estimating the lensing bending angle, we 
have assumed that we have a redshift survey that extends to a characteristic redshift z. Figure 
2 shows the results for large scale surveys of varying depths. In this figure, we have assumed 
the standard parameters for the MAP and SDSS projects. The characteristic depth of the SDSS 
redshift survey is z = 0.2. SDSS has photometric redshifts for 10® galaxies that should extend the 
survey to z = 1. These photometric redshifts are accurate to ±0.02 in redshift; certainly accurate 
enough to compute the projected surface densities needed to predict lensing. Figure 3 shows the 
cumulative signal-to-noise, where we have summed over all of the multipoles. We have divided 
(S/N) in equation (|26|) by a factor of 2 taking into account the limited sky coverage, rJsDSS- The 
predicted signal-to-noise is quite large (15 and 35) for the vacuum-dominated and standard CDM 
models. 

This cross-correlation, if detected, directly probes the gravitational potential fluctuations 
at low redshift. In principle, it should yield an accurate determination of the biasing factor. 
Armed with this measurement, we should be able to directly compare the gravitational potential 
fluctuations at decoupling with the gravitational potential fluctuations in the local universe. 

We thank Seljak and Zaldarriaga for providing their code to generate the intrinsic CMB power 
spectrum. DNS acknowledges the MAP/MIDEX project for support. MS and TS acknowledge 
support from Research Fellowships of the Japan Society for the Promotion of Science. 
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Fig. 2. — The cross-correlation, H{9), for ^fwhm = 0.21 degree. Solid and broken lines show the 
SCDM model and the ACDM model as in Fig. 1. For each model, the curves show the case in 

which Z = -Zdecj 

1, 0.5, and 0.2. 
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Fig. 3. — The cumulative signal-to-noise ratio. Solid and broken lines show the SCDM model and 
the ACDM model as in Fig. 1. For each model, the curves show the case in which z = z^ec, 1) 0-5, 
and 0.2. 
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